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Abstract. The local, uncorrelated multiplicative noises driving a second-order, purely noise-induced, ordering phase transi- 
tion (NIPT) were assumed to be Gaussian and white in the model of [Phys. Rev. Lett. 73, 3395 (1994)]. The potential scientific 
and technological interest of this phenomenon calls for a study of the effects of the noises' statistics and spectrum. This task 
is facilitated if these noises are dynamically generated by means of stochastic differential equations (SDE) driven by white 
noises. One such case is that of Ornstein-Uhlenbeck noises which are stationary, with Gaussian pdf and a variance reduced by 
the self-correlation time T, and whose effect on the NIPT phase diagram has been studied some time ago. Another such case is 
when the stationary pdf is a (colored) Tsallis' q-Gaussian which, being & fat-tail distribution for ^ > I and a compact-support 
one for q < \, allows for a controlled exploration of the effects of the departure from Gaussian statistics. As done before 
with stochastic resonance and other phenomena, we now exploit this tool to study — within a simple mean-field approximation 
and with an emphasis on the order parameter and the ''susceptibility' — the combined effect on NIPT of the noises' statistics 
and spectrum. Even for relatively small T, it is shown that whereas fat-tail noise distributions (^ > 1) counteract the effect of 
self-correlation, compact-support ones (^ < 1) enhance it. Also, an interesting effect on the susceptibility is seen in the last 
case. 
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INTRODUCTION 

As it is known, the role of temperature in equilibrium phase transitions can equally well be played by any set of 
spatially uncorrelated Gaussian white noises, regardless of their origin, provided that they have the same mean and 
variance and act additively on the system. Under the influence of multiplicative noises, the extended-system correlative 
of the phenomenon of noise-induced transitions (NIT) 1 1 1, namely a purely noise-induced phase transition (NIPT), may 
occur A comprehensive account of the many ways the phenomenon may take place can be found in |l2l|3]. However, 
for consistency with our previous work |l4','5l|, we shall restrict here to the 1994 model by Van den Broeck, Parrondo 
and Toral (VPT) |6 7 1, in which they proposed the following mechanism for the NIPT; 

1. An initially unimodal pdf gets rapidly destabilized towards a multimodal one. 

2. If spatial couphng is strong enough, the new states couple to form ordered domains that might subsequently 
coarsen and grow. 

Aiming at finding a nonequilibrium phase transition arising solely from the multiplicative nature of the noise, 
and characterized (in the limit of an infinite system) by ergodicity breakdown (only microstates compatible with the 
macroscopic broken symmetry should appear) and multiple steady state probability distributions P^^{{xi}), the authors 
in |6| set up a model that is the lattice version of a scalar reaction-diffusion model submitted to multiplicative local 
noises TJ,(f). The system's state at (continuous) time t is given by = U' stochastic variables Xi{t) defined at the sites 
r, of a hypercubic lattice, and obeying a system of coupled ordinary stochastic differential equations (SDE) 



f{xi)+g{xi)r]i + ^i, with f{x) = -x{l+x^f and g{x) = l+x^. 



(1) 



Here A,- = {D /2d)Y,jen{i){^j ~ ^i) is ^ discretization of the Laplacian V^;ic(r,r), D is the lattice version of the diffusion 
coefficient, and n(;) stands for the set of 2d nearest neighbors of r,. The Tj,(f) are uncorrelated, Gaussian and white 



(Tj,(f)> = , (ri,{t)iij{t^)) = a^dijd{t-t^). (2) 

As argued in H, establishing a phase transition rigorously is a difficult task even in equilibrium, where at least the 
explicit form of the steady-state pdf is known. Hence, one cannot resort to the traditional techniques from equilibrium 
statistical physics. However, the oldest and simplest ansatz that can reproduce (albeit not always faithfully) ergodicity 
breakdown, namely Weiss' mean field, can be readily adapted to this nonequilibrium situation. 

The method proposed by can be sketched as follows: By integrating the multivariable FPE over all variables 
except Xj and using the isotropy and translational invariance of the steady-state properties, one gets an exact (but 
implicit) steady-state equation for the one-site pdf in terms of the steady-state conditional average E{y) of y S 
n(/), given the value of x = x,. To determine the unknown function E{y), they introduce Weiss' mean-field (MF) 
approximation: to neglect fluctuations in neighboring sites, so that E{y) — x independent of y. The value of x then 
follows from the self-consistency condition 

x = m, with m = (x);f=const = J xP^^{x,x)dx = F{x). (3) 

When this nonUnear equation has multiple solutions, there are several corresponding steady state probabilities P^^{x) 
and the MF approximation predicts a phase transition with ergodicity breakdown (usually accompanied by symmetry 
breakdown). If, for example, / is odd and g even, then any realization {xi{t)} is equally probable as {— x,(f)} and one 
should expect {x) — 0. However, with the appearance of multiple solutions, this symmetry need not be fulfilled by the 
separate solutions, and one typically finds "ordered" phases with an order parameter m=\{x)\^ 0. 

Since F{x) is odd, i = is always a root of Eq. Hence, if the phase transition is second-order, its phase bound- 
ary — which provides rich qualitative information — lies where this root becomes unstable, i.e. where {dF jdx) \ = 1 . 
Now, the results of measurements and numerical simulations are correlations, order parameters and susceptibilities. 
Although the MF approximation is unable to predict the former ones, it yields predictions of the remaining two that 
can thus serve as a guide for numerical and (prospective) real experiments. The "susceptibility" we look upon is the 
MF correlative of the one defined in ||6l|71 for numerical simulations: 



1 



x^P^\x,x)dx — m"^ 



(4) 



Given the exploratory character of this work, we shall perform a still simpler MF approximation consisting in 
replacing in Eq. ([IJ A,- A,- = D{x — x, ), where i: is a parameter that will be determined self-consistently. Since the 
SDEs get decoupled, we hereafter consider a generic one: i = f{x) + g{x)ri + D{x — x) . We postpone the discussion 
of the numerical implementation of the MF method until we have described the process with q'-Gaussian distribution. 
For the sake of comparison with the results obtained in that case, we illustrate the findings in ||6l|7] by plotting in Fig. 
[l]the order parameter m and the "susceptibility" x as functions of the coupling D and the white noise intensity a . 
That the NIPT is a combined effect of noise and coupling is evidenced by the existence of a threshold value of D. Note 
also that the deterministic system is disordered however strong the coupling and that, of course, the NIPT is reentrant 
as a function of a^. 

Early in the long list of works inspired by 16] [71, two of us were involved in the study of the consequences of 
the multiplicative noise being Ornstein-Uhlenbeck (OU) ||4||3. In this case, Eq. ^ is replaced by (fj((f)) — and 
{'niit)lj{t')) = 5,;(c7V2T)exp(-|f-f'|/T), where the 77,(f) obey 

T77, = -77, + ^,-, with {^,{t)^j{t')) = a^5ij5{t~t'). (5) 

Although the colored noises {tj,} make the process {x,} non-Markovian, some approximations (interpolation schemes) 
render a Markovian (i.e. tractable) process, still capturing some of the essential features. Of course, at the price of 
adding a new unsystematic approximation to the MF one. However, we dispose of a neat control parameter, namely T, 
to compare with the white-noise case. Figure[2]is the corresponding plot for T = 0.1. For further details, see iHIS]. 
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FIGURE 1. Case of white multiplicative noise, a) order parameter and b) susceptibility as functions of coupling D and white-noise 
intensity . 




COMBINED EFFECT OF SPECTRUM AND STATISTICS 



One possible generalization of Eq. Q is 

d 



dri 



V,(77) + ^(f), with V,(t7) 



In 



1+ ^% ^ T? 



(6) 



proposed some time ago as model for coiTelated diffusion lH. As it occurred previously with the OU noise, this 
generalization provides a device to explore statistics effects by varying just one parameter (namely q, at constant T and 
ff^; note that the proper control parameter here is T{q — 1 ) /(J^, contrarily to the OU case in which it was t/cj^). 




a [arb. units] 



FIGURE 3. a) Detail of the unnormalized P^{t\) near its cutoff for q = 0.9, and b) phase boundaries in {a^,D) plane, for T = 0.1. 
Dotted line: q = 0.9; solid line: q= 1.0. dashed line: 151= 1.1. 



The stationary properties of the noise r\, including the time-correlation function, have been studied in ifTTIl so here 
we summarize the main results. The stationary probability distribution is given by 



1 



T 



^\ (7) 



where Zq is the normalization factor. This distribution can be normalized only for ^ < 3. The first moment (77) = is 
always equal to zero, and the second moment 



t(5-3^) 

2 

is finite only for <7 < 5/3, being larger than ^ for q>\. For q< \ the distribution has a cut-off, and it is only defined 
for |tj I < Tjf = ^ T-^'_q^ (Fig. 3 1). Finally, the correlation time Xq of the stationary regime of the process Tj (f ) diverges 

near q ~ 5/3 and it can be approximated over the whole range of values of q ss, XqK, 2t/(5 — 3>q). Clearly, when 
<7 ^ 1 we recover the limit of T] being a Gaussian colored noise, namely the Omstein-Uhlenbeck process ^of/(f), 

with correlations {^ou{t)^ou{t')) = fi^exp — |f — f'|/T and probability distribution ~ Z"' exp — ^^Qy. This 

process gives rise to interesting phenomena when it drives different kinds of nonlinear systems Ii9l [TOl [Tlf [121 [131 [T4l 

[I5][l6l? ,[i2i[ia. 

The expression for arises from a consistent Markovian approximation based on phase-space functional 

integration (details to be published elsewhere). The numerical implementation of the MF method compromises 
precision and speed: infinite integrals like those of Eqs. ([3| and Q [unless ^ < 1, see below] are performed by 
means of a 160 pt. Gauss-Hermite algorithm, whereas the finite integrals in the exponent of the stationary pdf are 
performed by means of a 96 pt. Gauss-Legendre algorithm. To determine the phase boundary, the Newton-Raphson 
algorithm is used (unless the corresponding function is badly conditioned, in which case a succession of finer sweeps 
is resorted to). In the ^ < 1 case, when the integrals over x in Eqs. Q and Q are bounded because of the bound in T], 
it would be faster to solve for the .x-bounds using the Newton-Raphson algorithm and then applying the 96 pt. Gauss- 
Legendre algorithm. However, the function seems not to be well conditioned for this algorithm. Hence we resort to 
naive integration inside a while loop. 

The integrals performed in the aforementioned way seem to be precise enough, except for large D/a^ where 
unphysical ordered states appear. For that reason, we have limited the exploration of the ((J^,D) plane to (a^ > 2) 
(except in Fig.[r| where the analytical expression for is used). Figures|4]and|5]exhibit the order parameter and 

susceptibility results respectively for ^ = 1.1 and q = 0.9 in the T = 0.1 case, showing that fat-tail noise distributions 
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FIGURE 4. Case of non-Gaussian multiplicative noise, with T = 0.1 and q = 1.1. a) order parameter and b) susceptibility as 
functions of coupling D and white-noise intensity CT^. 




FIGURE 5. Case of non-Gaussian multiplicative noise, with T = 0.1 and q = 0.9 a) order parameter and b) susceptibility as 
functions of coupling D and white-noise intensity CT^. 

> 1) counteract the effect self-correlation (namely, they advance the ordering boundary as is increased at 
constant D), and compact-support ones {q < 1) enhance it (they retard the ordering boundary). Particular interest 
rises the effect of {q < 1) multiplicative noises on the susceptibility: as seen in Fig.[5]3, it shifts from being larger on 
the ordering boundary to being larger on the disordering boundary. 

DISCUSSION 

As warned earlier, the character of this work is exploratory, and the value of its results is to be orientative of what to 
expect with numerical integration of the system of SDEs, as well as with more refined MF ansatze and/or interpolation 
schemes. In order not to run into the "forbidden" region, and given the ordering effect of ^ > 1 multiplicative noises, 
we have chosen to limit our exploration to « 10% around q= 1, just in order to discover trends. But in principle, our 
consistent Markovian approximation allows to explore the whole meaningful range of q values. 
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